Tiny stresses are capable of triggering earthquakes and tremors in Arunachal Himalaya

The Arunachal Himalaya has been hosting some notable events in the recent past. The tectonic history of Arunachal Himalaya is complex and has been influenced by several major tectonic events, including the 1950 Mw8.6 Assam-Tibet earthquake. In this study, we explored the effect of dynamic stresses generated by teleseismic events on the triggering of seismicity in the region. We analyzed 34 large teleseismic events since 2010 and found triggering during six events. The change in seismicity was also confirmed by analysis with the STA/LTA method. The triggering in the region occurred in the form of earthquakes and tremors. The dynamic stress as low as 1 kPa was found capable of triggering. The back-azimuth angle does not play an important role in the triggering. The angle direction of incoming waves with respect to the fault ~ 60° and ~ 120° is the possible reason for triggering in the region. The triggering occurred in the Mishmi and Main Central Thrust regions. The largest triggered event, ML2.3, was triggered 7.5 h after the 2012 Indian Ocean earthquake of Mw8.6. The region is tectonically very sensitive and tiny stresses are capable of triggering seismicity in Arunachal Pradesh.

www.nature.com/scientificreports/Himalaya, excluding Arunachal Pradesh.Hence, we are carrying out a systematic study of dynamic triggering in the Arunachal Himalayas.In this study, we attempted to identify the triggering in Arunachal Pradesh and explore the causes of the triggering.

Study region
Arunachal Himalaya covers an area between latitudes 26° 40′ N to 29° 25′ N and longitudes 91° 35′ to 97° 25′ E in the northeastern Himalayas.The northeastern Himalayas and its surrounding region are dominated by N-S oriented compressive tectonic stresses 28 .The region has two major's tectonic units viz. the Main Central Thrust (MCT) and the Main Boundary Thrust (MBT).The seismicity of the MBT bordering the lesser Himalayas and Tertiary layers is well-recognized, and numerous significant earthquakes have occurred along the boundary in the past 28 .The MCT is a significant intercontinental shear zone that separates the crystalline higher Himalayas in the north from the lesser Himalayas in the south.Other significant faults in the region comprise the Kopili Fault, Lohit Thrust, and Mishmi Thrus 29 (Fig. 1).Approximately 32 earthquakes of magnitudes ≥ 6 have occurred in 121 years (1900-2021) (ANSS database), which shows that the region is tectonically active and critically stressed.The region is also characterized by many hot springs and geysers 24 .One of the most prominent recent earthquakes was the M w 6.4 earthquake on January 4, 2005, with its epicenter near the India-Bhutan border.An event of M w 5.5 impacted the region near Assam-Arunachal Himalaya on April 24, 2019.The Arunachal Himalaya can be divided into the western and eastern regions by the Siang Window [30][31][32] .

Data and method
We selected 34 events that generated theoretical dynamic stresses ≥ 1 kPa at RUPA station.The theoretical dynamic stress in the area is computed 33 using Here T = 20 s and the value of A 20 are find out using the surface wave magnitude relation log 10 A 20 = M s − 1.66 log 10 Δ − 2, Where A 20 is in micrometers and Δ in degree 33 .The observed dynamic stresses are presented in Table 1.
The CSIR-NGRI, Hyderabad has been operating a seismic network of 34 broadband stations in the Arunachal Himalayas since 2010 (Fig. 1, Supplementary Table ST1).We examine the local catalog from years 2010 to 2018 comprising 848 events in Arunachal Pradesh.The local catalog's magnitude of completeness (M c ) was determined by using the maximum curvature method and found to be 2.3 (Fig. 2a).The M c is an essential element in figuring out data completeness 34,35 .The cumulative plot shows no important change in seismicity following all significant events (Fig. 2b).The catalog data typically lacks the smaller events that may be detected by waveform analysis; hence the waveform data give more information on the dynamic triggering than the catalog data 3,8,9,17,36 .
V ≈ 2πA 20 /T www.nature.com/scientificreports/After catalog analysis we analyzed the waveform data from this network to study the triggering.The waveform data of 24 h duration were analyzed to find the signature of triggering in the region.The instrument corrections were applied to the waveform data.The instrumental corrected components of N-S and E-W were rotated to a great circle path to obtain the transverse and radial components.The waveform data was high passed at 5 Hz to obtain the signature of local events since remotely triggered events generally have greater frequencies than the teleseismic event 2 .We generated the spectrogram of the waveform data to identify the triggered events using a short window Fourier transform.The waveform is analyzed manually to identify the triggered tremors or earthquakes.The statistical significance of the change in seismicity is computed using the β value.A longer duration of ~ 6 h of background seismicity is essential to estimate the β value, which can be estimated as 37,38 .
where N a and N are the number of events after the main shock in the targeted triggered window and the total number of events during the background and triggered window, respectively, T and T a are the total time window The events in Arunachal Pradesh that were examined for dynamic triggering.The events induced dynamic stresses of ≥ 1 kPa at the RUPA seismic station.For each station mean distance of the station from the event, range of backazimuth at all analyzed stations and the dynamic stress of each event were calculated.The regions peak dynamic stress may be calculated 14 using the peak ground velocity as σ = µ(PGV ) v , where σ, μ, PGV and v are dynamic stress, shear rigidity, Peak ground velocity, and surface wave phase velocity, respectively with an assumption of phase velocity of 3.5 km/s and shear rigidity of 35 GPa.www.nature.com/scientificreports/and duration of the triggered window respectively.A β value of 2 implies an increase in the seismicity rate, whereas a β value of -2 suggests a significant drop in the seismicity rate 13 .

Results
At the outset, during the examination of the waveform data of mainshocks, we used spectrogram and visual analysis to identify high-frequency microearthquakes and tremors during the large-amplitude surface waves of 34 teleseismic events and found triggering during six teleseismic events, (1) April The Banyak Islands, Sumatra earthquake occurred on April 6, 2010, at 22:15 UTC with a magnitude of M w 7.8 and thrust faulting mechanism at the plate boundary between the Australia-India and Sunda plates 39 .The peak dynamic stresses vary from 1.79 to 3.09 kPa with an average of 2.21 kPa.The regions peak dynamic stress calculated using the peak ground velocity as σ = (μ(PGV))/v, where σ, μ, PGV and v are dynamic stress, shear rigidity, Peak ground velocity, and surface wave phase velocity, respectively 14 with an assumption of phase velocity of 3.5 km/s and shear rigidity of 35 GPa.
The triggering is found at the ANNIG station which is 6 km away from the hot spring (Fig. 1) and peak dynamic stress was 3.09 kPa.The triggering phenomenon has been observed to occur within 300-400 s in the body wave phase and between 2000 and 3000 s in the surface wave passage (Fig. 3).The β value for the ANNIG station was 4.5 (Supplementary Fig. S1) which indicates a statistically significant increase in the seismicity.
The Mentawai earthquake on the western coast of Sumatra occurred at 14:42 UTC on October 25, 2010, M w 7.8 with a thrust faulting mechanism.It generated a large localized tsunami that hit the Mentawai Islands 40 .The Mentawai earthquake was recorded in 10 Indian stations and we find the triggering at RUPA station (Fig. 4).The triggered events are observed between 2200 and 3000 s during the surface waves and we recorded one microearthquake and three tremors during that duration.
The Tohoku earthquake occurred on March 11, 2011, M w 9.1 at 05:46:24 UTC, near the east coast of Honshu, Japan with a thrust faulting mechanism on the Pacific-North America sub-duction zone.The earthquake generated peak dynamic stresses of 12.4-16 kPa at the vertical component, averaging 13.97 kPa.The triggering is detected at the MIGN (Fig. 5) and RUPA (Fig. 6) stations.The triggering at the MIGN and RUPA stations was initiated during the surface waves.At the MIGN station, we observed the triggered tremors and at the RUPA station during the time period of 1700-3000 s, we observed one microearthquake and two tremors.Both the stations show evidence of delayed triggering (Figs. 5, 6).The β value at RUPA and MIGN station is 2.5 and 7.09, respectively.
The M w 8.6 earthquake in the eastern Indian Ocean on April 11, 2012, is the largest strike-slip event on record.This earthquake caused a sudden and significant increase in seismic activity worldwide 41 .The events were detected at the KAPU station in the form of instantaneous and delayed triggering (Fig. 7).We observed triggering in the form of microearthquakes and tremors.The KAPU station was also triggered during its largest aftershock of M w 8.2 (Fig. 8).The Indian Ocean earthquake and its largest aftershock generate β value 2.69  and 3.78 at the KAPU station.Specifically, two earthquakes with magnitudes of M L 1.8 and 2.3 were recorded at around 4.6 and 7.5 h, respectively after Indian Ocean main shock.The Iquique earthquake occurred on April 1, 2014, M w 8.2 at 23:46 UTC, off the coast of Chile.The earthquake occurred due to shallow depth thrust faulting between the Nazca and South American plates, where the Nazca plate is subducting at a rate of 65 mm/year 42 .The earthquake was recorded at 20 stations.The Iquique earthquake was triggered at stations JENG (Fig. 9), PANG, and KAPU (Supplementary Fig. S2, S3).At all 3 stations, we observed the triggering in the form of tremors.All the triggered stations are near the Mishmi thrust and Lohit thrust which are active faults 29,30 .The β value of triggered stations is > 2, which indicates a statistically significant increase in the seismicity.The remaining stations do not show any sign of triggering except the TAWG station which possibly showed triggered tremors (Supplementary Fig. S4).
The 2015 M w 7.9 Nepal earthquake and its aftershock of M w 7.3 generated dynamic stress of 201.5 and 180.6 kPa respectively.Despite having high dynamic stress, the Nepal earthquake was unable to trigger any seismicity in the study region (Fig. 10).Similarly, the earthquake with a magnitude of 6.8 occurred in Myanmar on August 24, 2016, with a peak dynamic stress of 62.6 kPa in the region, but unable to trigger events.The criticality of a region may be one of the criteria for triggering the seismicity in the region.
We do not find any triggering evidence on other stations for any events.We plotted the log envelope of nontriggering stations for the six events (Supplementary Fig. S6 to S10).We visually inspected all waveforms but could not find any sign of triggering at the rest of the stations.

Discussion
In the study, we systematically analyzed the dynamic triggering in the Arunachal Himalayas following large and shallow teleseismic earthquakes.Numerous studies provided evidence of dynamic triggering along major plate boundaries or volcanic/geothermal areas 8,14,43 and stable intraplate locations 19,22,23,44,45 .Examining dynamic triggering is a valuable tool for recognizing the appearance of earthquakes or tremors within a particular area due to minor stress perturbations.
Arunachal Himalaya is situated in a region prone to seismic activity owing to its proximity to the eastern segment of the Indian Plate boundary, where it comes into contact with the Eurasian Plate.This area is situated within the wider seismic zone of the Himalayan region and exhibits susceptibility to seismic activity.Earlier cases of dynamic triggering in the Himalayan region were found in south-central Tibet following the 2004 M w 9.1 Sumatra & 2005 M w 8.6 Nias earthquakes 26 , central Himalaya during the 2007 M w 8.5 Sumatra earthquake 27 , and the southwest China during the 2004 M w 9.1 Sumatra earthquake & 2012 M w 8.6 Indian Ocean earthquake 20 .
We examined 34 teleseismic events and identified the triggered event at six stations (ANNIG, RUPA, MIGN, KAPU, PANG, and JENG) following the six remote main events, namely, 2010 M w 7.8 Sumatra, 2010 M w 7.8 Mentawai, 2011 M w 9.1 Tohoku-Oki, 2012 M w 8.6 Indian Ocean, 2012 M w 8.2 Indian Ocean aftershock, and 2014 M w 8.2 Iquique earthquake.There are two basic hypotheses for understanding dynamic triggering: one based on the Coulomb failure criteria and the other linked with fluid movement or aseismic creep activation 13 .The delayed triggering is caused by fluid movement or aseismic creep, whereas the instantaneous triggering is caused by coulomb failure.Since we found both instantaneous and delayed triggering in our research, we believe that both models are responsible for triggering.It is interesting to note that previously a low-velocity zone (V s = 3.38 km/s) was found at a depth of 20-30 km by joint inversion indicating a weak mid-crustal zone 46 .Microearthquakes are generally found to occur at shallower depths ≤ 15 km in the region 47 .There are potential factors that may contribute to the occurrence of triggering in Arunachal Himalaya, with one of them being the presence of hot springs (Fig. 1). Figure 1 illustrates the proximity of ANNIG, TAWG, JENG, and PANG to the hot springs, which are seen as a potential factor in the initiation of triggering.Dynamic stresses cause fluids to migrate along fault pathways, altering the stress distribution within the Earth's crust; the fluid migration and the associated changes in stress can temporarily stabilize fault segments by increasing confining pressure and inhibiting immediate slip along the fault 48 .The triggered stations KAPU, JENG, and PANG are around 80-120 km and 20-80 km from the Mishmi and Lohit thrusts respectively (Fig. 1), both of which are seismically active 49 .The most recent strong earthquake in the Mishmi Thrust 50 was M w 6.5 in 2000.The KAPU station was also triggered following the 2012 Indian Ocean and its aftershock.
Many recent studies have shown that the direction of incoming waves could be a function of dynamic triggering 22,51,52 and have documented that the incidence angle of the triggering waves is parallel to the strike of the central ridge in Taiwan.Earlier studies identified that the faults in NE Iran parallel to the incoming waves are most likely to experience triggering 15 .Moreover, similar results have been found in the Coso Geothermal Region 53 .Our finding is consistent with the aforementioned observations.Hence, the incidence of the incoming surface waves relative to the faults is significant in the Arunachal Himalaya (Fig. 11).The events of Sumatra 2010, Mentawai 2010, and the Indian Ocean 2012 formed an angle of ~ 120° with regard to the fault located near the respective recording station.In a similar manner, the seismic events in Tohoku 2011 and Chile 2014 were shown to have a relative angle of ~ 60° with respect to the fault.Therefore, angles ~ 60° and ~ 120° play an important role in triggering (Fig. 11).
The back-azimuth angle (BAZ) is important in understanding how stress is conveyed.When seismic waves from a teleseismic event contact local geological features at a particular angle, they may cause stress to be transferred to neighbouring fault lines.The highest stress encountered during this interaction, represented by the peak dynamic stress, might increase the chance of local earthquakes being produced 15 .In our study, the triggering has been identified at various back azimuth angles (red circles) with respect to peak dynamic stress.We do not get a particular angle for the triggering (Supplementary Fig. S5).The highest dynamic stress was seen at an angle of 270°; however, no indication of triggering was observed at this specific angle.Therefore, it may be inferred that the back azimuth (BAZ) is an improbable factor for triggering in the region.a result of chemical interactions between different rock types and water 54 .The vibrations caused by wave energy act as a rapid jolt, stressing the cracks.These little cracks may start to spread because of the dynamic stress caused by seismic waves.Following the rules of subcritical fracture development, the cracks grow larger, which increases the stress at their points.The increased stress speeds up the pace at which these cracks expand, maybe to the point where they become large enough to cause a seismic event.Fundamentally, the interaction between dynamic forces and subcritical fracture formation creates a complicated process that may help to triggering in the region 54 .The dynamic stress threshold that triggers earthquakes varies according to geographical locations and the tectonic environment.Many recently conducted studies have shown that the dynamic stress of 1 kPa is sufficient to cause seismicity worldwide 55,56 .The dynamic stress range of 4-11 kPa can trigger the tremor in Japan.If we consider the case of the Himalayas we already discussed, there are only a few studies 20,26,27 related to dynamic triggering.In the central Himalayas, the ~ 9 kPa dynamic stress triggers approximately 40 earthquakes during the first 12 h of teleseismic wave arrival 27 .In the present study we found dynamic stress of 1 kPa is capable of triggering seismicity in this region.

Conclusion
A systematic investigation was conducted to examine dynamic triggering in Arunachal Himalaya using 34 significant remote and distant earthquakes between April 2010 and December 2018, having peak dynamic stress of at least 1 kPa.Using the waveform data, we identified triggered seismicity in the form of microearthquakes and non-volcanic tremors during six remote mainshocks.The triggering are not significantly influenced by the BAZ.The angle (~ 60° and ~ 120°) of incoming seismic waves with respect to the fault is a possible reason for the triggering.It is important to note that multiple processes may be occurring at the same time and contributing to the triggering in the region.The largest triggered event was M L 2.3, 7.5 h after the occurrence of the Indian Ocean earthquake, M w 8.6.The majority of the events with higher peak dynamic stress do not cause earthquakes or tremors.However, under some conditions, it is seen that a dynamic stress of 1 kPa is capable of triggering.We infer that the Arunachal Pradesh region is extremely stressed, and even tiny stresses can cause seismic triggering.

Figure 1 .
Figure 1.Tectonic map of Arunachal Himalaya and seismic stations (purple triangles), where green solid circle represents the hot spring.MCT: Main Central Thrust, MBT: Main Boundary Thrust, MFT: Main Frontal Thrust, BL: Bomdilla Lineament.The Seismotectonic data used for plotting the faults (solid orange line) and hot spring, Geological Survey of India, Government of India, Kolkata, India, last access on 09 October 2023, https:// bhuko sh.gsi.gov.in/ Bhuko sh/ Public.The figure is made using GMT version 6.3.0.

Figure 2 .Figure 3 .
Figure 2. (a) Plot of the log number of events shown in empty triangles and log of a cumulative number of events shown in empty circles versus magnitude.The red vertical line corresponds to the magnitude of the completeness of 2.3.(b) A temporal plot of the seismicity where variously coloured vertical lines depict the occurrence of trigger events.The cumulative number of events (solid orange line) does not indicate any substantial rise after the main shock arrival.

Figure 4 .
Figure 4. (a) Vertical components of the 5 Hz band passed waveform at RUPA station during the Mentawai earthquake, October 25, 2010, M w 7.8, (b), (c) and (d) are Zooming portion of blue vertical line in (a), (b) Transverse and Radial component, (c) Vertical component with a 5 Hz high pass filter, (d) Vertical component spectrogram, (e) zoom in the portion of the number marked on the spectrogram.Red and black colour waveforms represent (high pass 5 Hz) the triggered micro earthquakes and tremors in (e).(b) and (c) Follow the same time scale of (d).
https://doi.org/10.1038/s41598-023-49068-3www.nature.com/scientificreports/Another possible reason for triggering is subcritical crack growth.Dynamic stress is intimately related to subcritical crack growth, a progressive process where microscopic cracks inside a geological fault zone expand as

Figure 5 .
Figure 5. (a) Vertical components of the 5 Hz band passed waveform at MIGN station during the Tōhoku earthquake, March 11, 2011, M w 9.1, (b), (c) and (d) are Zooming portion of the blue vertical line in (a), (b) Transverse and Radial component, (c) Vertical component with a 5 Hz high pass filter, (d) Vertical component spectrogram, (e) zoom in the portion of the number marked on the spectrogram and black colour waveform represent (high pass 5 Hz) triggered tremors.(b) and (c) Follow the same time scale of (d).

Figure 6 .Figure 8 .
Figure 6.(a) Vertical components of the 5 Hz band passed waveform at RUPA station during the Tōhoku earthquake, March 11, 2011, M w 9.1, (b), (c) and (d) are Zooming portion of blue vertical line in (a), (b) Transverse and Radial component, (c) Vertical component with a 5 Hz high pass filter, (d) Vertical component spectrogram, (e) zoom in the portion of the number marked on the spectrogram.Red and black colour waveforms represent (high pass 5 Hz) the triggered microearthquakes and tremors in (e).(b) and (c) Follow the same time scale of (d).

Figure 9 .
Figure 9. (a) Vertical components of the 5 Hz band passed waveform at JENG station during the Iquique earthquake, April 1, 2014, M w 8.2, (b), (c) and (d) are Zooming portion of the blue vertical line in (a), (b) Transverse and Radial component, (c) Vertical component with a 5 Hz high pass filter, (d) Vertical component spectrogram, (e) zoom in the portion of the number marked on the spectrogram and black colour waveform represent (high pass 5 Hz) triggered tremors.(b) and (c) Follow the same time scale of (d).

Figure 10 .
Figure 10.(a) Vertical components of the 5 Hz band passed waveform at RUPA station during the Nepal earthquake, April 25, 2015, M w 7.8, (b) Transverse and Radial component, (c) Vertical component with a bandpass of 5 Hz, (d) Spectrogram of the vertical component.

Figure 11 .
Figure 11.Plot of incoming seismic waves at the triggered stations ANNIG, RUPA, KAPU, and JENG stations with respect to faults.The length of the arrows is proportional to the amplitude of the incoming surface waves, which is calculated as sqrt [(amplitude of Transverse) 2 + (amplitude of Radial) 2 ].The triggered and nontriggered teleseismic events (having dynamic stress ≥ 10 kPa) are shown by the red and yellow arrows.The Seismotectonic data used for plotting the faults (solid orange line), Geological Survey of India, Government of India, Kolkata, India, last access on 09 October 2023, https:// bhuko sh.gsi.gov.in/ Bhuko sh/ Public.The figure is made using GMT version 6.3.0.